PROCEEDINGS 

OF SCIENCE 



SFB/CPP-08-60, DESY 08-1 15, LTH 801, HU-EP-08/28 



OO 

O 

o 

(N 

< 

in 



Oh! 



> 
(N 

(N 

od 
O 
oo 
o 



Static-light meson masses from twisted mass lattice 
QCD 



Karl Jansen 

DESY, Platanenallee 6, D-15738 Zeuthen, Germany 



E-mail: k.arl.jansen@desy.de 



Chris Michael, Andrea Shindler 

Theoretical Physics Division, Department of Mathematical Sciences, University of Liverpool, 
Liverpool L69 3BX, UK 



E-mail: c . michael@liverpool .ac.uk 



E-mail: andrea . shindler@ liverpool .ac.uk 



Marc Wagner* 1 

Humboldt-Universitat zu Berlin, Institut fiir Physik, Newtonstrafie 15, D- 12489 Berlin, Germany 
E-mail: }ncwagner@physik . hu-berlin . de 




We compute the static-light meson spectrum using two-flavor Wilson twisted mass lattice 
QCD. We have considered five different values for the light quark mass corresponding to 
300MeV<mps <600MeV. We have extrapolated our results, to make predictions regarding the 
spectrum of B and B s mesons. 



The XXVI International Symposium on Lattice Field Theory 
July 14 - 19, 2008 
Williamsburg, Virginia, USA 



* Speaker. 

"'"On behalf of the ETM Collaboration. 



© Copyright owned by the author(s) under the terms of the Creative Commons AttribLition-NonCommercial-ShareAlike Licence. 



http: //pos . sissa . it / 



Static-light meson masses from twisted mass lattice QCD 



Marc Wagner 



1. Introduction 

A systematic way to study B and B s mesons from first principles is lattice QCD. Since amt > 1 
at currently available lattice spacings for large volume simulations, one way to proceed is to use for 
the b quark a formalism such as Heavy Quark Effective Theory (HQET). In this paper we consider 
the leading order of HQET, which is just the static limit. In this limit a "B meson" will be the 
"hydrogen atom" of QCD. States are either labeled by = (L± 1/2)^, where L denotes orbital 
angular momentum, ±1/2 the coupling of the light quark spin and & parity, or by S = (1/2)-, 
P = (1/2)+, P+ = (3/2)+, D = (3/2)", ... In the static limit these states will be doubly de- 
generate, since there is no interaction with the heavy quark spin, i.e. the total spin for given is 
F -'? =j-'?±\/2. 

The spectrum of static-light mesons has been studied comprehensively by lattice methods in 
the quenched approximation with a rather coarse lattice spacing [jl|]. More refined lattice studies 
with Nf = 2 flavors of dynamical sea quarks have also explored this spectrum [||, ||, |], ||, ^]. Here 
we also use Nf = 2 and are able to reach lighter sea-quark masses, which are closer to the physical 
ujd quark mass allowing more reliable extrapolations. 

2. Static-light meson creation operators 

In the continuum an operator creating a static-light meson with well defined quantum numbers 
is given by 



Q(x) I dnr^(n)U(x;x+m)xjf(x+rt). 



(2.1) 



Q creates an infinitely heavy antiquark, J dn denotes an integration over the unit sphere, U is a 
spatial parallel transporter, and y creates a light quark separated by a distance r from the antiquark. 

) is an appropriate combination of spherical harmonics and /-matrices yielding a well defined 
total angular momentum / and parity . The meson creation operators we use in this paper are 
listed in Table p} 





p3» 


r(^)(x) 


o h 


lattice J 3 " 


(1/2)- [5] 
(1/2)+ [P_] 


o-,i- 

0+,l+ 


1 , Yjxj 


Ax 
Ai 


(1/2)-, (7/2)-,... 
(1/2)+, (7/2)+,... 


(3/2)+ [P + ] 
(3/2)- 


l+,2+ 
l-,2" 


Yl%i — 72*2 (and cyclic) 
75(^1 - 72*2) (and cyclic) 


E 
E 


(3/2)+, (5/2)+,... 
(3/2)-, (5/2)-,... 


(5/2)- [D + ] 
(5/2)+ [F_] 


2-3- 
2+,3+ 


7l*2*3 + 72*3*1 + 73*1*2 
75 (7l*2*3 + 72*3*1 + 73*1*2) 


A 2 
A 2 


(5/2)-, (7/2)-,... 
(5/2)+, (7/2)+,... 



Table 1: static-light meson creation operators. 



When constructing lattice versions of the operators (2.1), one has to replace the integration 
over the unit sphere by a discrete sum over lattice sites with fixed distance r from the static anti- 
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quark. For J =1/2 and J = 3/2 we use six lattice sites, i.e. 



>(x) = (2(x) £ r (7 ^(n)C/(x;x + rn)yA(x + m), 
n=±ei ,±§2, ±63 



and for J = 5/2 we use eight lattice sites, i.e. 



<2(x) £ r (/ ^(n)£/(x;x + rn)i//tx + rn). 

n=±6i±e2±e3 



(2.2) 



(2.3) 



The spatial parallel transporters U in (2.2) are straight paths of links, while in (2.3) we use "diagonal 
links", which are averages over the six possible paths around a cube projected back to SU(3). The 
states created by these lattice meson creation operators do not form irreducible representations of 
the rotation group SO(3), but representations of its cubic subgroup Oh- Therefore, these states 
have no well defined total angular momentum /. They are linear superpositions of an infinite 
number of total angular momentum eigenstates. The common notation of the corresponding Oh 
representations together with their angular momentum content are also listed in Table |l|. Note that 
we do not consider Oh representations T\ and T 2 , because these representations yield correlation 
functions, which are numerically identical to those obtained from the operators listed. 



3. Simulation setup 

We use 24 3 x 48 gauge configurations produced by the European Twisted Mass Collaboration 
(ETMC). The fermion action is Nf = 2 Wilson twisted mass [0, |8[ and the gauge action is tree- 
level Symanzik improved JlO| ] with j3 = 3.9 corresponding to a lattice spacing a = 0.0855(5) fm. 
We consider five different values of the twisted mass /i q (cf. Table |J). Tuning to maximal twist has 
been performed at the lightest /i q value yielding K^ r = 0.160856 and "pion masses" in the range 
300MeV < m PS < 600 MeV. For details regarding this setup we refer to ^ ||]. 





raps in MeV 


number of gauge configurations 


0.0040 


314(2) 


1400 


0.0064 


391(1) 


1450 


0.0085 


448(1) 


1350 


0.0100 


485(1) 


350 


0.0150 


597(2) 


500 



Table 2: ji q values and corresponding pion masses mps. 



4. The static-light meson spectrum 

We determine the static-light meson spectrum from 6x6 correlation matrices. These matrices 
are built from trial states, which are generated by operators from Table [I] belonging to a fixed Oh 
representation. Due to parity breaking discretization errors of the Wilson twisted mass formulation 
we use states with = + and states with 8? = — in the same correlation matrix. 
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To reduce statistical noise and to simplify quark smearing, we use stochastic propagators ob- 



tained by inverting four random spin diluted timeslice sources per gauge configuration Q13|]. We 
use the HYP2 static action [jl4, 15] and apply Gaussian smearing [|^, |n|] to the dynamical quark 



operators with APE smeared spatial links [|18[]. For details regarding the construction and the com- 
putation of these correlation matrices we refer to @. 

Static-light meson masses are determined both by solving a generalized eigenvalue problem 
and computing effective masses, and by fitting a suitable ansatz of exponentials to the correla- 
tion matrices. Results obtained from both approaches are consistent. Note that static-light meson 
masses diverge in the continuum limit due to the self energy of the static quark. Therefore, we 
always consider mass differences of static light mesons, where this self energy exactly cancels, and 
which are physically meaningful in the continuum limit. Mass differences between various states 
and the lightest static-light meson (J^ = (l/2)~ ground state) are shown in Figure [I] as functions 
of (mps) 2 . 
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Figure 1: static-light mass differences linearly extrapolated to the u/d quark mass and the s quark mass. 



We extrapolate linearly in (mps) 2 to the physical u/d quark mass (nips = 135 MeV) as well 
as to the physical s quark mass (taken here as mps = 700 MeV). Note that we consider the unitary 
sector in both cases, where valence quarks and sea quarks are of the same mass. This implies for 
the s quark extrapolation a sea of two degenerate s instead of a sea of u and d. We plan to improve 
these computations by using ATy = 2+1 + 1 flavor gauge configurations currently produced by 
ETMC [20]. Results are shown in Figure [T] and Table || together with the corresponding # 2 /dof 
indicating that a straight line is a suitable ansatz. 

Performing a similar u/d extrapolation for the P wave mass difference yields 
m{P + ) = 140(22)MeV, i.e. the (1/2)+ state is lighter than the (3/2)+ state as is usually 

expected. We see evidence for a reversal of this level ordering as the light-quark mass increases ob- 
taining a difference m{P-) —m{P + ) = 98(33) MeV at the strange quark mass. It will be interesting 
to check this in the continuum limit. 
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u/d quark extrapolation: 

/ t ^ \ / r-t\ ■ » i it 

m(J^) - m(S) m MeV 


s quark extrapolation: 
m(J^) - m(S) in MeV 


lit c 

X /dof 


(1/2)-* 


[S*] 


678(40) 


994(90) 


1.01 


(1/2)+ 


[P-] 


354(23) 


574(38) 


0.20 


(3/2)+ 


[P+] 


495(11) 


473(20) 


2.05 


(3/2)" 


[D-] 


807(27) 


885(48) 


0.04 


(5/2)" 


[D + ] 


876(28) 


874(55) 


2.81 


(5/2)+ 


[F-] 


1154(39) 


1232(67) 


0.72 



Table 3: static-light mass differences linearly extrapolated to the u/d quark mass and the s quark mass. 



5. Predictions for B and B s meson masses 

To make predictions regarding the spectrum of B and B s mesons, we interpolate between our 



static4ight lattice results and experimental data for charmed mesons pip. To this end we assume 
a linear dependence in 1/rriQ, where thq is the mass of the heavy quark. Results for P wave B and 
B s states are shown in Figure || and Table ^. Note that the lines labeled by "S(J"^ = (l/2)~)" in 
Figure || together with the experimental values for B* and B* (the triangles intersected by these 
lines) indicate that straight lines are suited for interpolation. 
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Figure 2: static-light mass differences linearly interpolated to a heavy b quark, a) B mesons, b) B s mesons. 
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In Table |] we compare our lattice results to experimental data available for P wave B and 

|j. The lattice results for B states are larger by around 15% than 



B s mesons [21, 22, 23, 24 



their experimental counterparts, while there is significantly better agreement for B s mesons. For a 
conclusive comparison it will be necessary to investigate the continuum limit, which amounts to 
considering other values for the lattice spacing. 
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m(F& = 


0-)in MeV 




m - 


-m(F& = 


0") in MeV 


state 


lattice 
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C0 


state 


lattice 


PDG* 


CDF 


C0 


B 


400(22) 








B s o 


507(27) 








B* 


416(22) 


T 








572(27) 


T 






B\ 


513(8) 


419(8) 


455(5) 


441(4) 




478(14) 


486(16) 


463(1) 




B* 2 


524(8) 


I 


459(6) 


468(4) 


b; 2 


512(14) 


I 


473(3) 


473(3) 



Table 4: lattice and experimental results for P wave B and B s states (*: unknown). 



In Figure |b we also plot the BK and B*K threshold (406 MeV and 452 MeV respectively). Our 
P wave B s results indicate that the corresponding decays are energetically allowed. Consequently, 
one can expect that these states have a rather large width compared to e.g. certain excited D s 
mesons. 



6. Conclusions 

We have studied the static-light meson spectrum by means of Nf = 2 Wilson twisted mass 
lattice QCD. We have considered five different values of the dynamical quark mass corresponding 
to 300MeV<mps < 600 MeV. We have performed an extrapolation in the light quark mass to the 
physical u/d mass and s mass respectively, as well as an interpolation in the heavy quark mass to 
the physical b mass. Our results agree within <15% with currently available experimental results 
for P wave B and B s states. 

Future plans regarding this project include an investigation of the continuum limit, which 
amounts to considering other values for the lattice spacing. Moreover, we plan to perform similar 
computations on Nf = 2 + 1 + 1 flavor gauge configurations currently produced by ETMC. In 
particular, this will allow us to include a sea of u and d in B s computations. We also intend to 
determine the static-light decay constants fs and fs s . 
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